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ABSTRACT 

Aims. We present a spectroscopic and photometric analysis of the rapid rotator KIC 8429280, discovered by ourselves as a very young 
star and observed by the NASA Kepler mission, designed to determine its activity level, spot distribution, and differential rotation. 
Methods. We use ground-based data, such as high-resolution spectroscopy and multicolor broad-band photometry, to derive stellar 
parameters (v sin i, spectral type, T^, log g, and [Fe/H]), and we adopt a spectral subtraction technique to highlight the strong chro- 
mospheric emission in the cores of hydrogen Ha and Ca n H&K and infrared triplet (IRT) lines. We then fit a robust spot model to the 
high-precision Kepler photometry spanning 138 days. Model selection and parameter estimation is performed in a Bayesian manner 
using a Markov chain Monte Carlo method. 

Results. We find that KIC 8429280 is a cool (K2 V) star with an age of about 50 Myr, based on its lithium content, that has passed 
its T Tau phase and is spinning up approaching the ZAMS on its radiative track. Its high level of chromospheric activity is clearly 
indicated by the strong radiative losses in Can H&K and IRT, Ho-, and H/3 lines. Furthermore, its Balmer decrement and the flux 
ratio of Can IRT lines imply that these lines are mainly formed in optically-thick sources analogue to solar plages. The analysis 
of the Kepler data uncovers evidence of at least seven enduring spots. Since the star's inclination is rather high - nearly 70° - the 
assignment of the spots to either the northern or southern hemisphere is not unambiguous. We find at least three solutions with nearly 
the same level of residuals. Even in the case of seven spots, the fit is far from being perfect. Owing to the exceptional precision of the 
Kepler photometry, it is not possible to reach the noise floor without strongly enhancing the degrees of freedom and, consequently, 
' the non-uniqueness of the solution. The distribution of the active regions is such that the spots are located around three latitude belts, 

i.e. around the star's equator and around ±(50° -60°), with the high-latitude spots rotating slower than the low-latitude ones. The 
equator-to-pole differential rotation d£l = 0.27 rad d _1 is at variance with some recent mean-field models of differential rotation in 
rapidly rotating main-sequence stars, which predict a much smaller latitudinal shear. Our results are consistent with the scenario of a 
higher differential rotation, which changes along the magnetic cycle, as proposed by other models. 

> ' 

OO ' Key words. Stars: activity - stars: starspots - stars: rotation - stars: chromospheres - stars: individual: KIC 8429280 - X-rays: stars 

^J. , 1 . Introduction low latitudes, often of opposite polarity in the two hemispheres 

\£) • (Nels on et al.i 2011). On the other hand, theoretical models of 
. Starspots are characteristics of solar-like activity observed in ration of magnet ic fields in fully convective pre-main se- 

; cool stars They are tracers of magnetic flux tube emergence and e (pMS) starfj afgue ^ ±e bafjic d Q action ifj fa ^ 

. can provide information about the different forces acting on the ff2 ime ( Kuker & Rmzei 1999). 
■ • ■ flux tubes during their buoyant rise, and on photospheric mo- 
. > • tions, such as the latitudinal drift of spots during the activity The differential rotation pattern is expected to be dictated by 

V ! cycle and the differential rotation, which are both essential in- the Taylor-Proudman theorem, so that the iso-contour lines are 

5_i ■ gredients of mean field dynamo theory. cylinder-shaped. 

However, even though the activity cycle in Sun-like stars Recent calculations of a 1-Mq rapidly rotating main- 
is believed to be dr iven by the meridio nal circulation with a sequence star with period P = 1.33 days, predict a value of 0.08 
low eddy diffusivity dBonanno et al.ll2002|) . the situation is much radd 1 , surprisingly close to the Solar value despit e there being 
less clear for late-type rapidly rotating stars. In this latter case, a fact or of 20 between the average rotation rates (Kuke r et al.l 
numerical simulations suggest that the magnetic dynamo ac- 1201 ll) . However, for HD 171488, a young Sun with an equato- 
tion can produce wreaths of strong toroidal magnetic field at rial rotation period of 1 .33 days, values of differe ntial rotation as 

high as 0.50 rad d" 1 ha ve been found (e.g., iMarsden et al.ll2.Q06l: 
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Based on public Kepler data, on observations made with the Italian 



A unique opportunity for trying to solve these open problems 



Telescopio Nazionale Galileo operated on the island of La Palma by is 8 iven b ? the NAS A Ke P ler mission ' which is Providing an un- 

the Centro Galileo Galilei of INAF (Istituto Nazionale di Astrofisica) at precedented data set of the photometric variability of a large Star 

the Spanish Observatorio del Roque del los Muchachos of the Instituto sample (|Borucki et al.||2010|). Its combination of very high pho- 

de Astrofisica de Canarias, and on observations collected at the Catania tometric precision and long, uninterrupted time coverage is es- 

Astrophysical Observatory (Italy). sential not only for its primary goal (the discovery of transiting 
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exoplanets) but also for the study of many astrophysical phe- 
nomena in stellar variability, such as the magnetic activity and 
differential rotation in cool stars. 

Single active stars in the field can be most efficiently selected 
on the basis of their high coronal emission, and turn out to be 
mostly young stars with an age of a few hundred Myr, i.e. in 
the zero-age main sequence (ZAMS), or even younger, namely 
weak T Tauri (wTT S) or post-T Tauri (pTTS) stars (see, e.g., 
iGuillout et al.l l2009.and reference therein). 

KIC 8429280 (=TYC 3146-35-1) was selected as an X- 
ray active star from the cross-correlat i on of the ROS AT All- 
Sky Survey (RASS; IVoees et all 1 19991 l2000l) with Tycho and 
Hipparcos catalogs dPerrvman et all Il997l) . It is a relatively 
bright star (V = 9.9 mag) and falls in the field of view of Kepler. 



2. Ground-based observations and data reduction 

2.1. Spectroscopy 

Three spectra of KIC 8429280 were acquired with SARG, the 
echelle spectrograph at the Italian TNG telescope (La Palma, 
Spain). The first one was taken on 2007 May 31 within a sur- 
vey of optical counterparts to X-ray sources (proposal TAC67- 
AOT15/07A). With the yellow grism and a slit width of a/8, this 
spectrum has a resolution of R = A/AA 57 000, covering the 
spectral range 4600-7900 A in 55 echelle orders, with a signal- 
to-noise ratio S/N~70. The other two spectra (S/N ranging from 
about 60 to 100) were taken, adopting the same slit width, on 
2009 Aug 11 (proposal TAC71-AOT20/09B) with the red and 
blue grisms and cover the wavelength ranges 5500-1 1 000 A and 
3600-5 100 A, respectively. 

Another spectrum of KIC 8429280 with S/N-50 was taken 
at the M.G.Fracastoro station (Serra La Nave, Mt. Etna, 1750 
m a.s.l.) of the Osservatorio Astrofisico di Catania (OAC - Italy) 
on 2009 July 13. The 91 -cm telescope of the OAC was equipped 
with FRESCO, a fiber-fed echelle spectrograph that covers the 
spectral range 4300-6800 A in 20 orders with a resolution R 
21000. 

For calibration purposes, spectra of radial and rotational ve- 
locity standard stars (see Table[T]i, as well as bias, flat-field, and 
arc-lamp exposures were acquired at both observatories during 
each observing run. 

The data reduction was performed with the echelle task 
of the IRAJ-Q packa ge, following standard steps (see, e.g., 
ICatanzaro etalll2010l) . 

2.2. Photometry 

CCD images in the Johnson-Cousins B, V, Rc, and Iq filters were 
acquired using the focal -reducer imaging camera with the 91 -cm 
telescope of the OAC on 2010 November 1 1 . Data reduction was 
carried out following standard steps of overscan region subtrac- 
tion, master-bias subtraction, and division by average twilight 
flat-field images. The BVRqIc magnitudes were extracted from 
the corrected images through aperture photometry performed 
with DAOPHOT using the IDLE! routine APER. Standard stars 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of the Universities 
for Research in Astronomy, inc. (AURA) under cooperative agreement 
with the National Science Foundation. 

2 IDL (Interactive Data Language) is a registered trademark of ITT 
Visual Information Solutions. 



Table 1. Radial/rotational velocity standard stars. 



Name 


Sp. Type 


A V 


v sin i d 


iNoies 






(Km s ; 


( Km s ) 




HD 115404 


K2V 


7.60" 


3.3 


RV 


HD 145675 


K0V 




0.8 


v sin i 


HD 157214 


GOV 


-79.2 b 


1.6 


v sin i 


HD 10700 


G8V 


-17. l b 


0.9 


v sin i 


HD 221354 


Kl V 


-25.20 3 


0.6 


RV, vsin; 


HD 32923 


G4V 


20.50 11 


1.5 


R V, v sin i 


HD 12929 


K2III 


-14.6 C 


1.6 


RV 


HD 187691 


F8V 


-0.2 C 


2.8 


RV 


HD 182572 


G8 IV 


-100.35 11 


1.9 


RV, vsin; 


HD 161096 


K2III 


-12.53 3 


2.1 


RV 



a l Udrv et all dl999h. b lNordstr6m et all d2004h. 
c lEvansUl967h . " iGlebocki & Gnacinskil l2005). 



in the cluster NGC 7790 (IStetsonll2000l) were observed to calcu- 
late the transformation coefficients to the Johnson-Cousins sys- 
tem. The zero points for the V magnitude and B - V, V - R, and 
R — I colors were determined by the obser vation of a standar d 
star close to KIC 8429280 (BD+42 3339. ICastelaz et ai1ll99ll) . 
Photometric data are summarized in Table [2] 




6160 6170 6180 6190 



Fig. 1. Continuum-normalized observed spectrum (solid black 
line) in the region between 5220 and 5240 A {top panel) and 
between 6100 and 6200 A (middle and bottom panels) of 
KIC 8429280. The synthetic spectrum calculated with the pa- 
rameters listed in Table[3]is superimposed (smooth red line). 
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Table 2. Johnson-Cousins photometric data. Uncertainties are in parenthesis. 



Name 


B V 


Rc 


Ic 


r 


H-' 




KIC 8429280 


10.88 (0.04) 9.94 (0.04) 


9.38 (0.05) 


8.86 (0.06) 


8.114 (0.021) 


7.672 (0.031) 


7.541 (0.018) 



From 2MASS catalog jCutri et all2003h . 



3. Target characterization from ground-based 
observations 

3. 1 . Astrophysical parameters 

High-resolution spectroscopy permits us to derive most of 
the fundamental stellar parameters, such as radial (RV) and 
projected rotational velocity (vsim), spectral type, luminosity 
class, effective temperature (T e s), gravity (logg), and metallic- 
ity([Fe/H]). 

Cross-correlation functions (CCFs) were computed with the 
IRAF task fxcor to derive both RV and vsim. For this pur- 
pose, as templates we used spectra of slowly rotating RV stan- 
dard stars listed in TableQ] We average d the results from ind ivid- 
ual echelle orders as described, e.g., in Fra sca et"a"D ([2010). The 
CCFs were used for the determination of vsim through a cali- 
bration of the full-width at half maximum (FWHM) of the CCF 
peak as a function of the v sin i of artificially broadened spec- 
tra of slowly rotating standard stars (TableQ]) acquired with the 
same setu p and in the same o bserving seasons as KIC 8429280 
(see, e.g JGuillout et alJl2009h . obtaining a value of 39.3 ±2.4 
km s _1 , where the error is the standard deviation of v sin i values 
from each individual echelle order. 

We measured, within the errors, the same RV in the SARG 
spectra of KIC 8429280 taken in 2007 and 2009 (RV = 
-33.1 ± 0.5kms _1 , and -33.1 ± 0.4 km s -1 , respectively) and 
the FRESCO spectrum of 2009 (RV = -32.2 ± 0.8 km s 1 ). The 
constant RV measured by ourselves and the absence of any sign 
of duplicity in both the CCF and the line spectrum, suggest that 
KIC 8429280 could be a single star, although more spectra and 
high angular resolution images are needed to safely exclude the 
presence of a stellar companion (see Sect. 15.11 for a wider dis- 
cussion). 

We used the ROTFIT code dFrasca et al.l2003l l2006h to eval- 
uate r e ff, logg, [Fe/H], and re-determine vsin/. We adopted, 
as reference spectra, a library of 185 ELODIE Archive spec- 
tra of standard stars well distributed in effective temperature, 
spectral t ype, and gravity, and in a suitable range of metal- 
licities dPrugniel & Soubiranl 1200 lb . The SARG spectrum at 
Rsarg - 57 000 taken with the yellow grism was convolved with 

a Gaussian kernel of width W = A ^jl /^ EL0DIE - 1 /^sarg ^ t0 
match the resolution of ELODIE spectra (^elodie = 42 000). 
We applied the ROTFIT code to the echelle orders with a fairly 
good S/N, which cover the range 4600— 6800 A, in the observed 
spectrum. For each order, we estimated the stellar parameters of 
KIC 8429280 as the mean of the parameters of the ten reference 
stars that most closely resemble the target spectrum, which is 
quantified by a x 2 measure, and adopted their standard devia- 
tion as a measure of the uncertainty. The adopted estimates for 
the stellar parameters come from a weighted mean of the val- 
ues for all the individual orders, the weight accounting also for 
the x 1 (more weight to the most closely fitted or higher S/N or- 
ders) and the "amount of information" contained in each spectral 
region expressed by the total spectral-line absorption. The stan- 



dard error in the weighted mean was adopted as the uncertainty 
estimate for the final values of the stellar parameters. 

For the v sin i determination, we ran ROTFIT using as tem- 
plates a smaller grid of SARG spectra of inactive and slowly 
rotating standard stars (Table[TJl acquired with the same instru- 
mental setup as for our target. This allowed us to remove any 
eventual systematic error in the artificial broadening of template 
spectra introduced by the different resolution and provided us 
with a measure of v sin i independent of the CCF method. The 
vsin/ determined with ROTFIT, 37.1 ±2.9 kms _I , agrees with 
the value derived from the CCF. For these particular tasks, the 
Ha and Li i lines, as well as spectral regions affected by telluric 
features, were excluded from the analyzed spectral range. 

We also determined the stellar parameters by fitting selected 
spectral regions with sy nthetic spectra ge nerated by SY NTHE 
dKurucz & Avrettl 119811) using ATLAS9 dKuruczl 1 1 993h atmo- 
sphere models. When applying all the models, we considered 
a solar opacity distribution function and microturbulence veloc- 
ity f = 1 km s _1 , which are typical of a star with T e g = 5000 K 
and log g-4.5 (see, e.g.. lAllende Prieto et 1012004 . 

As starting values of T e ff and logg, we used those derived 
from ROTFIT. At the same time, we determined the projected ro- 
tational velocity by matching the metal lines present in our spec- 
tral range. We a dopted the lists of spectra l lines and the atomic 
parameters from lCastelli & H ubrig (2004), who updated the pa- 
rameters listed orig inally by Kurucz & Bell ( 1995). The stellar 
parameters were determined by minimizing the difference be- 
tween the observed and the synthetic spectra, using the x 1 as the 
goodness-of-fit parameter. Errors were estimated to be the vari- 
ation in the parameters that increases the^ 2 of a unit. 

In Fig. [TJ we show the result of our fitting of two different 
wavelength regions containing important neutral iron and cal- 
cium lines, which are well suited to constraining temperature, 
gravity, and metallicity. The cores of the two strong Fe i lines at 
/15227.2 and /15232.9 are possibly influenced by NLTE effects 
and chromospheric filling. 

In Table [3] we list the stellar parameters derived with both 
ROTFIT and SYNTHE, which agree very well with each other. 

3.2. Chromospheric activity and lithium abundance 

We evaluated the level of chromospheric activity from the emis- 
sion in the cores of the Can H&K and IRT lines, as well 
as in the Balmer Ha, Hp, and He lines. The star's age was 
inferred from the lithium abundance. We estimated both the 
lithium equivalent width and the chromospheric emission level 
with the "spectral subtrac tion" technique (see, e.g.. lHerb"igT l985: 
iFrasca & Catalano1ll994l) . This technique is based on the sub- 
traction of a "non-active", "lithium-poor" template made with 
observed spectra of slowly rotating stars of the same spectral 
type as KIC 8429280 with a negligible level of chromospheric 
activity and no apparent lithium absorption line. The equivalent 
width of the lithium line (Wu) and the net equivalent width of the 
Ca n H&K, Ca n IRT, Ha, and H/3 lines (W em ) were measured, in 
the spectrum obtained after subtracting from the target the non- 
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Table 3. Astrophysical parameters. Uncertainties are in parenthesis. 



Name Sp. Type 


RV 




l°g£ 


[Fe/H] 


v sin 


W u 


logJV(Li) 




(kms- 1 ) 


(K) 






(kms- 1 ) 


(mA) 




KIC 8429280 K2V 


-33.1 (0.5) 


5055 (135) a 


4.41 (0.25) a 


-0.02 (0.10)" 


37 (3) a 


270 (20) 


2.9 (0.1) 






5000(100) b 


4.50 (0.10) b 


-0.05 (0.10) b 


39 (3) b 







a From ROTFIT. 
b FromSYNTHE. 



Call H & K 




3928 3930 3932 3934 3936 3938 3940 
A (A) 



3962 3964 3966 3968 3970 3972 3974 
A (A) 



Fig. 2. Top panel: Continuum-normalized S ARG spectrum of KIC 8429280 (solid line) in the Ca n H & K region together with the 
inactive stellar template (dotted red line) broadened at the vsin i of the target and Doppler-shifted according to the RV difference. 
Bottom panels: The difference spectrum. He emission, already apparent in the observed spectrum, is clearly emphasized in the 
residuals. 



active template, by integrating the residual emission (absorption 
for the lithium line) profile. 

The Ca ii H & K lines display strong emission cores that ap- 
pear nearly symmetric without any evidence of self-absorption. 
The He emission is clearly visible and becomes more evident 
after the subtraction of the non-active template (Fig. [2). The 
peak intensity of the H&K lines is similar to t hat displayed 
by SA O 51891, a very young K0-1 star studied bv lBiazzo et alJ 
d2009h . 

In Fig. [3] the spectral subtraction reveals Ca n IRT lines with 
emission reversals in their cores and a remarkable filling-in of 
the HB line. Furthermore, the Ha line is completely filled-in by 
emission, with high values of net equivalent width (see TableHJ). 
This behavior is typical of very active stars and the Ha pro- 
file is quite similar to those frequently displayed by LQ Hya, 



a well-studied young star fo r which the Ha changes from filled- 
in to weak emission profi les dFekel et al.ll98 6: Strassmeie r et alJ 
ll993UFrasca et al.ll2008l) . The similarity between KIC 8429280 
and LQ Hya goes beyond the Ha profile, since the spectral 
type is the same for both stars (K2 V), the effective temperature 
is nearly identical, and the rotation period of LQ Hya is only 
slightly longer (1.6 versus 1.2 days). 

We also evaluated the total radiativ e losses in the chr omo- 
spheric lines following the guidelines of lFrasca et alJ d2010l) . i.e. 
by multiplying the average W em by the continuum surface flux at 
the wavelength of the line. The latter was evaluated by means of 
the spectrophotometric atlas of iGunn & Str vker ( 198 3J) and the 
angul ar diameters calculated by applying the Barnes & Evans! 
( 1 19761) relation. The net equivalent widths and the chromo- 
spheric line fluxes are reported in Table [4] 
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Ha Call XX8498.8542 




W V 




JL 





0.5 



0.6 



0.2 



0.0 



4850 4855 4860 4865 4870 4875 
A (A) 



6550 



6560 6570 
A (A) 



6580 



8500 8520 8540 
A (A) 



8560 



Call A8662 




8650 8660 8670 8680 
A (A) 



Fig. 3. Top of each panel: Observed, continuum-normalized SARG spectra of KIC 8429280 (solid line) in the HB, Ha, and Can 
IRT regions together with the non-active stellar template (dotted red line). Bottom of each panel: Difference between observed and 
template spectra. The residual Ha profile is plotted shifted downwards by 0.2 for the sake of clarity. The hatched areas represent the 
excess emissions that have been integrated to get the net equivalent widths. 



Table 4. Line equivalent widths and radiative chromospheric 
losses. 



Line 


Date 




Error 


Flux 




(yyyy/mm/dd) 


(A) 


(A) 


(ergcirT 2 s -1 ) 


Ho- 


2007/05/31 


0.922 


0.107 


3.35x10 s 




2009/07/13 


1.026 


0.150 


3.73x10 s 




2009/08/11 


0.999 


0.071 


3.63x10 s 


ti/' 


2009/07/13 


0.315 


0.118 


1.16x10 s 




2009/08/11 


0.442 


0.076 


1.63x10 s 


He 




0.264 


0.040 


0.48x10 s 


CanH 




1.075 


0.068 


1.66x10 s 


CanK 




0.795 


0.045 


1.45x10 s 


Can IRT ^8498 




0.538 


0.065 


1.45x10 s 


Can IRT A8542 




0.624 


0.065 


1.69x10 s 


Can IRT ,18662 




0.596 


0.074 


1.63x10 s 



On the basis of the Ha and HB flux, we evaluated a Balmer 
decrement Fua/Fu/s in the range 2.2-3.2. Values of the Balmer 
decrement in the range 1-2 are typical of optically thick emis- 
sion by solar and stellar plages (e.g., Chester] 1199 It iBuzasil 
Il989h ." while prominences seen off-limb give rise to values of 
~ 10. The flux ratio of two Can-IRT lines, ^8542/^8498 = 1.2, 
indicative of high optic al depths, is in the range of the values 
found by Chester ( 1991) in solar plages. Solar prominences have 
instead values of ^8542/^8498 ~ 9, typical of an optically-thin 
emission source. 

This suggests that the bulk of chromospheric emission of 
KIC 8429280, in both the Balmer and the Can lines, is basi- 
cally due to magnetic surface regions similar to solar plages, and 
eventually prominences will play a marginal role. 

The subtraction method also permits the measurement of 
the lithium equivalent width cleaned up from the contamina- 
tion of the close Fei A 6707.4 A photospheric line (Fig. @). 
The lithium equivalent width, Wu = 270 + 20 mA, abov e 
the Pleiades (lOOMyr) upper envelope dSoderblom et al.llT993h . 
translates into a very high lithium abunda nce, log N(Li) g 2.9 + 
0.1, ad opting the calibrations proposed bv lPavlenko & Magazzul 
(Il996h . With this lithium abundance and T eS = 5055 K, the 
star falls on the upper envelope of the a Per cluster (50Myr, 
ISestito & Randichll2 005). indicating an age of * 50Myr. 



For comparison. iFekel et al.l (Q986) measured for LQ Hya a 
slightly lower lithium equivalent width, Wu = 234 mA, deduc- 
ing a lithium abundance log N(Li) ^2.8, and concluded that that 
star should be at least as young as the Pleiades. This supports the 
measurement of a very young age for KIC 8429280, which could 
still be in the wTTS or pTTS phase and not yet have reached the 
ZAMS. 



0.6 















V 








_[ 1 


Co 1 i 











6700 



6710 6720 
A (A) 



6730 



Fig. 4. Continuum-normalized spectrum (solid black line) in the 
lithium region of KIC 8429280 together with the template spec- 
trum (dotted red line) of the non-active lithium-poor reference 
star, broadened at the v sin i of the target and Doppler-shifted ac- 
cording to the RV difference. 



3.3. Spectral energy distribution 

The standard BVRqIc photometry (Tabled, c omplemented wit h 
JHK S magnitudes from the 2MASS catalog (ICutri et al.ll2003l) . 
allowed us to reconstruct the spectral energy distribution (SED) 
from the optical to the near infrared (IR) domain. 

To perform a fit to the SED, we used the grid of Ne xtGen 
low-resolution synthetic spectra (Ha uschildt et al.l fl999) with 
logg = 4.0 and 4.5 and solar metallicity. The effective tem- 
perature was set to 5000 K in agreement with the value found 
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KIC 8429280 

Photosphere model 

□ Synthetic fluxes weighted with the BVRJJHK, bonds 
• BVRJJHK, dereddened fluxes 



till' :'t9 



~ -8.6 - ^ 



ii 



V. 



v 



4.0 
log(A [A]) 



E 0.010 



Frequency (doys ) 



Fig. 5. Spectral energy distribution for KIC 8429280 (dots). The 
NextGen spectrum with T e g = 5000 K and log g = 4.0 is over- 
plotted with a continuous line. 

with ROTFIT and SYNTHE (Tabled. The interstellar extinction 
(Ay) was evaluated from the distance according to the rat e of 0.8 
mag/kpc found by Miko laiewska & Mikolaiewskil d!980t) for the 
sky region around CI Cyg (very close to our target). Finally, the 
angular diameter (<p), which scales the model sur face flux over 
the ste llar flux at Earth, was allowed to vary. The ICardelli et al.l 
( 1989) extinction law with Ry = 3.1 was used. The best-fit solu- 
tion was found by minimizing the x 2 of the fit to the BVRqIcJ 
data, which are dominated by the photospheric flux of the star 
and are normally not appreciably affected by infrared excesses. 
The angular diameter derived from the SED, <p = 0.1 40 mas, im- 
plies a distance d ^ 50 pc if we adopt a ZAMS radius of 0.75 Rq. 
Adopting a radius of 0.88 Rq, which is typical of a 50-Myr old 
star with the same temperature as KIC 8429280 (see Sect. l4.U . 
a distance of about 60 pc is derived. 

As seen in Fig. [5] the SED is quite well reproduced by the 
synthetic spectrum and no excess is visible at near-IR wave- 
lengths. 

This provides a further check of the effective temperature 
and assures that this star has by far exceeded the T Tau phase, 
largely dissipating its accretion disk. However, a thinner "debris 
disk", whose presence could be only detected by mid- or far-IR 
data (lacking for KIC 8429280), cannot be excluded. 

4. Kepler light curve and spot modeling 

4.1. Photometric data 

We analyzed the long-cadence data of quarter (Q0) to quarter 
2 (Q2). The 6174 data points are taken from 2009 May 2 to 2009 
September 16, spanning nearly 138 days with a sampling of one 
photometric point each 30 minutes. The most pronounced gap, 
between Ql and Q2, amounts to about 4.5 days. 

The power spectrum of the Kepler time series, displayed in 
Fig. [6] shows two main peaks very close in frequency (0.828 
and 0.863 d~'), which correspond to periods of about 1.207 and 
1.159 days, respectively. The period uncertainty, estimated from 
the FWHM of the spectral window, is about 0.006 days. The 
two low-amplitude peaks at frequency of as 1.66 and 1.72 d 
are overtones of the two main peaks. 

The double-peaked periodogram is a clear fingerpr int of dif- 
ferential rotation, as shown, e.g., bv lLanza et al. I (fl99l in their 
simulation of light curves of differentially rotating spotted stars. 



Fig. 6. Power spectrum of the Kepler Q0+Q1+Q2 time series of 
KIC 8429280. 

They show that a Fourier analysis of uninterrupted time series 
of photometry with high precision (^f - 10~ 5 - 10~ 6 , which is 
comparable with the precision of the Kepler data) may detect a 
Sun-like latitudinal differential rotation. 

A parameter needed to apply a spot model to our data is the 
inclination of the rotation axis with respect to the line of sight. 
From the values of v sin i, stellar radius R, and rotation period P 
the inclination of the rotation axis can be estimated as 



The stellar radius cannot be estimated from effective temper- 
ature and luminosity, because this star has no Hipparcos par- 
allax and the Tycho value is unreliable. From Eq. [TJ a lower 
limit follows, 0.85 Rq, which exceeds the radius of a ZAMS 
star with T eB = 5055K and a mass of 0.76 M©, R a 0.75 fl©. 
Thus, we considered a radius of a younger star at the same tem- 
peratu re. We found, from the evolutionary tracks of Siess et al. 
(2000), that a ~50Myr star of 0.9 Mq has the same temperature 
of our target, but a radius of 0.88 Rq. With this value, we find 
sin; ^ 0.965, which implies an inclination of about 75°. Such 
a high inclination is confirmed by the value of i ^ 70° derived 
from the light-curve solution (see Sect. 14. 3t . 

4.2. Bayesian photometric imaging 

Photometric imaging can be a many parameter problem. An effi- 
cient and robust method is light-curve fitting using circular spots. 
We chose to use Dorren's (1987) analytical starspot model in its 
generalization to a quadratic limb-darkening law. Any departure 
from the spherical shape caused by the star's fast rotation was 
neglected. 

In the case of seven spots, we have 46 free parameters. Three 
parameters describe the star: inclination angle ;, equatorial angu- 
lar velocity Q eq , and equator-to-pole differential rotation d£l. To 
reduce the number of free parameters and enforce a monotonic 
dependence on latitude, (3, differential rotation is prescribed by a 
Sun-like sin 2 -law: 

Q(y8) = Q eq -t/Qsin 2 /?. (2) 

An additional parameter describes the spot rest intensity k. For 
each spot, six parameters are needed: the longitudes in an arbi- 
trarily rotating reference system at the beginning and at the end 
of the observation, the spot area at a given time, and three further 
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parameters describing spot area evolution (for details see below). 
The tw o coefficients u + an d u~ in the quadratic limb-darkening 
law (cf . lAlonso et al.ll2.Q08l) were fixed using the tabulate d coef - 
ficients from stellar model atmosphere calculations (Sing 2010): 
u + = 0.7042 and uT = 0.3774. The limb darkening is assumed 
to be the same for the unperturbed photosphere as well as for the 
spots. 

The parameterized law of differential rotation (Eq. [2]) pro- 
vides the absolute value of a spot's latitude /3 from its rotational 
frequency. To find out the hemisphere to which a spot belongs, 
we made a test run with 128 Markov chains in parallel, trying out 
all possibilities to arrange seven spots in the two hemispheres. 
Since the inclination was allowed to take negative values, there 
are always two Markov chains describing the same spot allo- 
cation situation. From this 64-hour test run, the three different 
solutions mentioned below (Sect. l4.3l l emerged. It cannot be ex- 
cluded for certain that there are further solutions. 

Spot evolution is described by a simple ansatz: The spot 
area, expressed in units of the star's cross-section, is assumed 
to evolve basically linearly with time. For a decaying spot, this 
prescription allows in principle an estimate of the turbulent mag- 
netic diffusivity. The spot evolution model even allows for one 
sudden change in the slope of the area-versus-time relation. 
Hence, altogether four parameters are needed to describe the 
area of a spot and its evolution: area at the time of the (sharp) 
bend, this time itself and two time derivatives of the spot area. 

In addition to the free parameters of the model, the derived 
parameters, whose physical meaning is relevant, are considered. 
We are interested in their marginal distributions, in order to de- 
rive the mean values and error bars. Examples are the latitude 
and period of a spot. The latter simply follows from the initial 
and final longitude of the spot, the former (except for its sign, 
i.e. the stellar hemisphere on which the spot is located) from the 
spot's rotation period and the two coefficients of the sin 2 -law for 
the differential rotation. 

To estimate all these parameters from the data, a Bayesian 
method was applied. This provides a straightforward way of es- 
timating parameters and their uncertainties from the data alone. 
Before analyzing the data, proper prior distributions have to be 
assigned. If information about the star's inclination i is lack- 
ing and selection effects are ignored, cos i should be evenly dis- 
tributed for < i < n/2. All non-dimensionless parameters such 
as radii or periods are represented by their logarithms. This en- 
sures that the posterior distribution for a radius will be consis- 
tent with that of an area and likewise the posterior for a period 
with that of a frequency, i. e. it makes no difference whether one 
prefers radii or areas, periods or frequencies. 

The likelihood function is constructed as follows. Spotted 
stars that are geometrically similar exhibit the same light curve, 
except for an offset in magnitude. This offset is considered ar- 
bitrary and removed by integration. Only the shape of the light 
curve matters. Hence, it is unnecessary to specify the magni- 
tude of the unspotted star, which is anyhow unknown. It is also 
assumed that the measurement errors have a Gaussian distribu- 
tion. The integration over the offset can then be done analyti- 
cally, which is quite convenient. 

With the N data points di gained at times t\, their standard de- 
viation cr, the fit fo(ti,Pi ■ ■ -Pm) and an offset co, the likelihood 
reads 



n 



A(o-,c (h pi . ■ .pM',di) = 

(di-fo(ti,pi . ..p M ) - c ) 2 



■ exp 



2cr 2 



The M unknown parameters are denoted by pj, with j — 1 . . . M. 
Integrating analytically away first the measurement error cr - us- 
ing Jeffreys' l/<x-prior (cf. Kass & Wassermanl ll996h - and 
then the uninteresting offset co, one gets a likelihood depend- 
ing on spot modeling parameters p\ . . . Pm only. It is this mean 
likelihood, averaged appropriately over measurement error and 
offset, from which the posterior density distributions for all un- 
knowns are obtained by marginalization. We can argue that the 
method itself determines error and offset from the data. 

When the light curve consists of different parts, first the like- 
lihood of each individual part is determined. Afterwards all like- 
lihoods are multiplied. This means, each part has been assigned 
its own error level and offset. 

We note that, besides an offset, a linear trend in the data may 
be taken into consideration, where the slope is also integrated 
away analytically. 

The Markov chain Monte Carlo (MCMC) method (cf. 
iPress et ail 120071) has been applied to explore the likelihood 
mountain in the high-dimensional parameter space. Basically, 
MCMC is able to find the globally best solution. 

4.3. Results 

There are six gaps exceeding a few hours in the combined light 
curve. Accordingly, the light curve has been divided into seven 
parts, with each part being assigned its own offset and error level. 
Hence, artificial vertical jumps between the data sets do not mat- 
ter. 

The MCMC method found three solutions, all with the same 
level of residuals. Two of the three solutions are quickly relaxed. 
The non-relaxed solution needs in any case very dark spots with 
spot intensity around 0.2 and would lead to a rather high value 
of the equator-to-pole differential rotation of 0.35 radd -1 . We, 
therefore, dismiss that solution. Nevertheless, all three solutions 
show nearly the same spot periods, which means that the max- 
imal observed rotation shear is remarkably the same: 0.21 
radd 1 . What is different are the latitudes of the spots and the 
spot evolution. In what follows, we prefer the solution with the 
lowest degree of differential rotation. All figures refer to that so- 
lution if not stated otherwise. Table provides values of some 
parameters for both relaxed solutions. 

A quick glance at the light curve (Fig. |7]i reveals the pres- 
ence of two enduring principal periods. The rotational periods 
are indeed tightly grouped around P 1.16 and P =s 1.20 days 
(Fig. [8]i as already indicated by the power spectrum (Fig. [6). 

The stellar surface obviously displays a plethora of long last- 
ing features. Within our framework of up to seven dark circular 
spots, we are unable to fit the light curve down to the measure- 
ment errors. At the most, we can fit a seven-spot model to the 
data, with the residuals being 2.36 mmag rms. The residuals, ex- 
ceeding by far the measurement errors, are undoubtedly due to 
the incapability of the model to cope fully with the data. 

The gain in goodness-of-fit relative to, e.g., a six-spot model 
is striking. To give an example, in a first analysis the combined 
Q0 and Ql data sets only had been considered. In the case of 
this 45-day subset, the residuals shrink from 1.26 mmag to 1.04 
mmag if one enlarges the number of spots from six to seven. 
The Bayesian information criterion (BIC) clearly jus tifies the 
increa se in the number of free parameters. The BIC dSchwarzl 
1 19781) allows theory ranking taking into account the goodness- 
of-fit as well as the number of free parameters used. With N, the 
number of measurements, cr 2 , the mean variance in one measure- 
ment, and M, the number of free parameters, the criterion reads 
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Fig. 7. Robust fit of the Kepler light curves. Q0 is from 953.53825 to 963.24461, Ql from 964.51 154 to 997.98296, and Q2 from 
1002.76455 to 1091.46678. The solid line is the best fit of the first solution. The residuals, shown at the top, are ±2.36 mmag. The 
gaps are artificially bridged by straight connecting lines. 



KIC8429280: spot periods 
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Fig. 8. Marginal distributions of spot periods. The two slow- 
est revolving spots exhibit the same period. The periods group 
around two principal periods, which produces the obvious "beat- 
ing" phenomenon quite understandable. 



BIC = N ■ log(cr 2 ) + M ■ \og{N). It expresses Occam's razor in 
mathematical terms. A model with a low BIC should preferably 
be used. 

Attempts to reduce the number of free parameters without 
any loss of information (according to BIC) were in vain. For ex- 
ample, we tried to fit the light curve with two different periods 
only. In all cases, the loss of goodness-of-fit, as expressed by o~, 
outclasses the gain in credibility by far. Thus, one has to accept 
that - owing to the unprecedented precision of the data - it is 
impossible to reduce the number of free parameters noticeably. 
That all these free parameters are well constrained by the data 
is testified by the mere fact that the Markov chains generally 
quickly converge to a relaxed solution. On the other hand, we 
did not try to reduce the size of the residuals of the fit by en- 
hancing the complexity of the spot model (more spots, different 
spot temperatures, faculae, etc.), because this would require an 
increase in the number of free parameters with a strong growth 
of the computing time and a loss of solution uniqueness. 

The marginal distribution for the inclination is shown in 
Fig. [9] The inclination angle is photometrically very well con- 



KIC8429280: inclination 




i [dig] 

Fig. 9. The star's inclination from Kepler photometry. Vertical 
lines mark the mean and the 68-per-cent confidence region. The 
inclination value confirms the relative large radius of the star. 

strained, as is apparent from the smallness of the 68-per-cent 
confidence region. 

The spot intensity, measured in units of the unspotted photo- 
sphere, is =! 0.57 (cf. Fig. [Toll which sounds reasonable, because, 
in the case of su nspots, the bolometric contrast is 0.67 (e.g., 
lLanza et al. 2003,not distinguishing between umbra and penum- 
bra). For KIC 8429280, a bolometric contrast of 0.57 would 
translate into a temperature difference of 650 K. 

The spots are rather large and they slowly evolve as appears 
in Fig. [TJJ Despite their large areas, the spots are never really 
touching each other. Hence, spot evolution does not seem to be 
influenced by the need to avoid spot overlapping. One should 
note that the spot evolution ansatz makes sense within the time 
span of photometric observation only. It must not to be extrapo- 
lated to any later times. 

As already mentioned, in terms of spot area evolution, the 
two solutions are quite different. We, therefore, provide for the 
second solution the corresponding spot area evolution plot, too 

(Fig.ua. 

A major difference between the two solutions concerns spots 
5 and 7. Spot 7 of the first solution (the largest one at the end) 
has partly replaced spot 5 of the second solution. This change 
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KIC8429280: spot intensity 
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Fig. 10. Spot intensity with respect to the unspotted photo- 
sphere. Vertical lines mark the mean and the 68-per-cent con- 
fidence region. 
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Fig. 11. Evolution of spot areas. Areas are expressed in units 
of the star's cross-section. The largest spots are those at high 
latitudes. 



in identity is already indicated in the periods. One should also 
note that there are, apart from spot 3, small but nevertheless 
highly significant differences in the periods of all the other spots 
(Tablegjl. 

That there are at least two solutions is not surprising. The 
inclination is rather high, so it is difficult to determine to which 
hemisphere a spot belongs. 

The main result concerns differential rotation. In the case of 
KIC 8429280, the spots are obviously long lasting, allowing us 
to separate the effects of spot evolution from differential rota- 
tion. However, even then a cautionary note is in place because 
of the large number of spots involved. In principle, the MCMC 
method should be able to find all possible spot periods, in prac- 
tice one gets at best one feasible set of periods - a conceivable 
"scenario", not more! The assignment of some dip in the light 
curve to a certain spot is not unambiguous. We lack the informa- 
tion about what happens when a spot is out of sight. 

The marginal distribution of the equator-to-pole differen- 
tial rotation is given in Fig. [13] For our preferred model, it is 
0.266 + 0.001 radd 1 , given a sin 2 law as a basis for the latitu- 
dinal dependence of differential surface rotation. The second so- 



KIC8429280: spot area evolution 
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Fig. 12. Evolution of spot areas for the second solution. The 
different types of line have the same meaning as in Fig.fTTl 

KIC8429280: differential rotation 
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Fig. 13. The equator-to-pole differential rotation. Vertical lines 
mark the mean and the 68-per-cent confidence region. 



lution results in an even larger differential rotation of 0.297 ^'002 
radd" 1 . 

As an additional check of the differential rotation, we tried to 
fit the Kepler light curve with a different set of free parameters, 
i.e. by allowing the spot latitudes to vary without any constraint 
on the differential rotation law. To save computing time, we ap- 
plied this model to the Q0+Q1 dataset, which is, however, large 
enough to constrain the rotation period of each spot. As a result, 
we found that the slowly revolving spots are at high latitudes. 
Despite the spot latitudes being notoriously ill-defined by pho- 
tometry, at least for ground-based data, the stellar surface rotates 
beyond any doubt in a Sun-like way. 

Expectation values with l-cr confidence limits for some in- 
teresting parameters are given in Table|5] The reader should be 
aware that the estimated parameter values and their (often sur- 
prisingly small) error bars are those of the model constrained by 
the data. The error bars simply indicate the "elbow room" of the 
model, nothing more. 

There is a pronounced long-term trend in the light curve. 
The star gets dimmer with time. To prevent any instrumental 
effect affecting the outcome of the photometric spot modeling, 
we allowed not only for offsets, but, moreover, for adjustable 
slopes in the seven time chunks. These slopes are considered 
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Table 5. Two 7-spot solutions: listed are expectation values and 
1-cr uncertainties. Latitudes f3 are derived from the assumed law 
of differential rotation (Eq. |2j. Periods P are given in days, the 
spot intensity k is in units of the intensity of the unspotted sur- 
face. The differential rotation c/Q (rad d~') is the equator-to-pole 
value of the shear. 



parameter 


1st solution 


2nd solution 


inclination 


i 




69!5 


+0.4 
-0.4 




69!6 


+0.4 
-0.4 


1st 


latitude 






62!5 


+0.2 
-0.2 




5TA 


+0.5 
-0.3 


2nd 


latitude 






-0?5 


+0.4 
-0.2 




-0?8 


+0.6 
-0.2 


3rd 


latitude 






62?5 


+0.2 
-0.2 




57!5 


+0.5 
-0.3 


4th 


latitude 


A 




17!3 


+0.2 
-0.2 




16?5 


+0.2 
-0.2 


5th 


latitude 






-53!9 


+0.2 
-0.2 




56!7 


+0.5 
-0.3 


6th 


latitude 






-8?0 


+0.2 
-0.2 




10!4 


+0.3 
-0.3 


7th 


latitude 


07 




-58!9 


+0.2 
-0.2 




-49! 1 


+0.2 
-0.4 


1st 


period 


Pi 


1 


20430 


+0.00004 
-0.00004 


1 


20390 


+0.00004 
-0.00004 


2nd 


period 


Pi 


1 


15793 


+0.00007 
-0.00007 


1 


15741 


+0.00009 
-0.00009 


3rd 


period 


Pi 


1 


20430 


+0.00005 
-0.00005 


1 


20430 


+0.00006 
-0.00007 


4th 


period 


P 4 


1 


16298 


+0.00007 
-0.00007 


1 


16253 


+0.00006 
-0.00005 


5th 


period 


Ps 


1 


19617 


+0.00012 
-0.00012 


1 


20345 


+0.00004 
-0.00004 


6th 


period 


Pi 


1 


15903 


+0.00005 
-0.00005 


1 


15948 


+0.00003 
-0.00003 


7th 


period 


Pi 


1 


20107 


+0.00008 
-0.00007 


1 


19479 


+0.00013 
-0.00014 


spot intensity 


K 




0.574 


+0.006 
-0.006 




0.599 


+0.017 
-0.015 


diff. rotation 


dQ 


0.2656 


+0.0011 
-0.0012 


0.2974 


+0.0016 
-0.0032 



The skewed error estimates for the latitude of the second spot are 
calculated under the assumption that it is located south of the equa- 
tor. 



nuisance parameters and analytically integrated away. The fit is, 
of course, tighter, with the residuals now being reduced to 2.17 
mag rms. Most importantly, our main results are unaffected. The 
inclination is lowered by merely three degrees, the spot intensity 
is slightly increased to 0.635, and the equator-to-pole differential 
rotation is increased somewhat to 0.269 + 0.003 radd 1 , which is 
between the values of the two solutions set out in Table[5] Owing 
to the smallness of these changes, there is no need to be doubtful 
about the long-term stability of the data. 

An additional improvement could be achieved by including 
more free parameters in the spot area evolution ansatz. The easi- 
est way to do this is to substitute the sharp bend between the two 
slopes with a more gradual transition. This can be achieved, for 
example, by a linear interpolation between the values of the two 
slopes of the first derivative of the area-versus-time relation over 
some time span, which is the new free parameter. The longer 
this time, the smoother the bend between the two time evolu- 
tion regimes. With seven free parameters more, the residuals are 
reduced from 2.36 mmag to 2.25 mmag, which is, because of 
the large number of data points involved, formally significant, 
but the main outcome remains nevertheless unchanged. The in- 
clination has decreased by a half of degree to i = 68°9 + 0.4 
and the equator-to-pole differential rotation has decreased to 
0.254 + 0.001 radd 1 . 



5. Discussion 

5.1. Is KIC 8429280 a binary? 

The constant radial velocity and the absence of secondary peaks 
in the CCF (Sect. [3Tj seem to exclude that KIC 8429280 is a 
spectroscopic binary in a short-perio d orbit. 

However, in the WDS Catalog dWorlev & Doudassl [19971) 
the source is classified as a close visual binary composed of two 
stars of equal magnitude with a separation of 073 and a position 
angle PA = 22°. The two components are resolved only in one 
visual obs ervation perfor med in 1988 with the 0.5-m Nice re- 
fractor by Cou teaul(ll990l) . Howev er, although the star is present 
in the Tycho Double Star Catalog (Fabricius et al. 2002), no en- 
try for the separation and position angle is present, even though 
in this catalog binaries with separations p > 072 and compo- 
nents of similar magnitude are normally resolved. Evidently, the 
Tycho instrument on-board the Hipparcos satellite, which was 
operational from December 1989 to February 1993 (only 1-5 
years later than the Couteau observation), was unable to distin- 
guish the components of this possible binary. As far as we know, 
no other hint of binarity is present in the literature. 

If we consider two identical stars with the same luminos- 
ity, temperature (T ef{ = 5055 K), and radius (R = 0.88 Rq), the 
distance estimated from the SED increases to about 85 pc. At 
this distance, an apparent separation of 073 corresponds to about 
25 AU. From the third Kepler law, with a mass of about 0.9 Mq 
per each component, we can estimate a minimum orbital period 
^orb *93 years, assuming that the semi-major axis of the or- 
bit is a a sin i = 25 AU. Thus, given the long orbital period, 
the apparent separation of the two components should have been 
nearly the same during the Couteaus's and Tycho observations. 

Moreover, in this hypothesis we can evaluate the semi- 
amplitude of the RV curve as 



k = 



2n 



P orh (l-e 2 Y' 2 (l-e 2 ) 1 ' 2 



kms 



(3) 



were we assume that a sin i = 25 AU. An RV variation should be 
easily detectable with observations taken a few years apart. 

From all the previous considerations, it is unlikely that 
KIC 8429280 is a double star. However, more spectra and high 
angular resolution images are needed to safely exclude the pres- 
ence of a stellar companion. 

5.2. Inclination 

The photometrically obtained value of the inclination, i a 69.5°, 
and the spectroscopically determined v sin i allow one to deter- 
mine the radius of the star. Assuming the shortest rotational pe- 
riod to be associated with the equator, one arrives at R — 0.93 R G , 
which is close to the expected value (see Sect. 13. 31 . As the incli- 
nation depends crucially on the model assumptions, this corre- 
spondence strengthens our confidence in the reliability of our 
spot-modeling approach. 



5.3. Spot contrast 

In a first analysis, considering quarters Q0 and Ql only, we 
identified rather large spots with a very low spot-to-photosphere 
contrast {k a 0.8). After analyzing a three-times longer time 
span, the spot surface brightness (Fig. ITDb in the Kepler pho- 
tometric band amounts to the comforting value of 60 % of the 
photospheric intensity. The corresponding temperature contrast 
r s p/7p h is 0.87. 
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In the case of 4 PMS stars in Orion (~2 Myr old), 
iFrasca et al. I d2009) found a temperature contrast Tsp/Tph in the 
range 0.7 1-0.92. For fairl y young stars (~ 20 0^1-00 Myr) such as 
a: 1 Cet, HD 166, and e Eri. Biazzo et all (120071) derived values of 
r s p/r p h ~0.85 from the simultaneous modeling of temperature 
(from line-depth ratio analysis) and light curves. The spots of 
KIC8429280 are, therefore, comparable to both those of Orion 
PMS and mildly active stars. However, we recall that a "spot" 
can be more than a single dark spot, namely an active region, 
which is a mixture of dark and bright surface features. 

However, despite their very high precision, the Kepler mag- 
nitudes are taken in "white light", i.e. with only one very broad 
band. Thus, there is no possibility of constraining the spot 
contrast without a contemporaneous multi-band photometric or 
spectroscopic monitoring. 

5.4. Differential rotation 

One of the unexpected results of the present work is the degree 
of differential rotation of KIC 8429280. Assuming a sin 2 law for 
the latitudinal dependence, one arrives at 0.27-0.30 radd for 
the equator-to-pole difference in angular velocity. This is about 
fiv e times larger t han in the case of the Sun. 

iBarnes et al. I J200l showed that in general differential ro- 
tation increases with stellar effective temperature. The highest 
value reported by th em for stars a s cool as ~ 5000 K, dQ » 0.2, 
is that measur e d by [Donati et alj d2003l) for LQ Hya in 2000. 
iMarsden et al.l d201 1 ) extended this dataset with values of dQ 
(in the range 0.08-0.45 radd -1 ) measured only for stars more 
massive than the Sun. 

For LQHya, which has similar stellar parameters and 
is rotating a little bit slower than our target (cf. Sect l3.2l ). 
iKovari et all {2004) determined the differential surface rotation 
from the cross-correlation of several Doppler maps. A solar-type 
differential rotation law, i.e. the equator rotating faster than the 
poles, with dQ. = 0.022 radd 1 (l ap time of » 280 day s) is found 
by the aforementioned authors. Dona ti et"ai] ((2003) measured 
surface differential rotation in LQ Hya and AB Dor by means of 
a new Doppler imaging technique that includes among its free 
parameters those of a solar-type differential rotation law. Their 
results show that the equatorial angular velocity and the ampli- 
tude of the latitudinal shear of AB Dor and LQ Hya change on 
timescales of a few years. This is ascribed to the redistribution 
of the angular momentum inside the stars likely to be caused by 
the action of a non-linear hydromagnetic dynamo. In particular, 
they found for LQHya a decrease in the shear between years 
2000 (dQ. 0.20 radd 1 ) and 2001 (0.01 radd 1 ). 

HD 141943, another very young more-massive star that is 
still in the PMS phase, has values of dQ ranging f rom about 0.23 
to 0.44 radd -1 in different epochs as derived by IMarsden et alj 
(1201 lh . 

Very different values of differential rotation have also been 
found for HD 171488, a young (-50 Myr) Sun. For this star, a 
very high solar-type differential rotation dQ « 0.4-0.5 radd -1 , 
with the equator lapping t he pol es eve ry 12-16 days, w as found 
by both IMarsden et all d2006l) and Ueffers & DonatH d2008l) . 
Much weak er values (dQ » 0.04 ) we re derived instead fo r the 
same star bv lJarvinen et al .] d2008h and lKovari etafl d2010h . 
A weak differential rotati on (dQ between 0.017 and 0.056 



rad d 1 ) was found for e Eri by Frohlichlj 2007 ) from a Bayesian 
reanalysis of the MOST light curve icrolll l2006t ICroll et alJ 



2006).' This is a mildly active star with the same spectral type as 
KIC 8429280 (K2 V), but with a considerably longer rotational 
period (P lot =s 11.2 days). 



Thus, from these observations, it seems that a higher dif- 
ferential rotation is more frequently encountered in rapid rota- 
tors, although different indications come from previous works, 
mainly based on ground-based photometry, which seem to dis- 
play the opposite behavior, i. e. a differential rotation d ecreasing 
with the rotation period (e.g jMessina & Guinanl200"l . The pre- 
cision of the ground-based light curves does not allow to draw 
firm conclusions and accurate photometry from space, as well as 
Doppler imaging, is strongly needed to enlarge the database of 
differential rotation measurements that must be compared with 
the results of new dynamo models for stars in different evolu- 
tionary phases (including the PMS stage) and with both low and 
high rotation rates. 

The high differential rotation that we found for KIC 8429280 
disagr ees with the predictions of the model of iKiiker et"al] 
d201 lh . who instead found rather low values of dQ w 0.08 for 
an (evolved) solar-mass star rotating with a period as short as 
1.3 days. However, this specific model was made for a star with 
a remarkably different mass, effective temperature, and inner 
str ucture. 

ICovas et"aT] d2005l) developed dynamo models capable of re- 
producing variations in surface differential rotation along the ac- 
tivity cycle, although they cannot produce anything as extreme 
as the episodes of almost rigid surface rotation (dQ * 0.0) as 
reported for LQ Hya. On the other hand, these models can re- 
produce variations in the differential rotation of several percent 
for stars with deep conve ctive envelopes rotating ten times faster 
than the Sun. lLanzal d2006l) found that the large value of the rota- 
tion shear (~ 0.2) observed for LQ Hya in the year 2000 implies 
a dissipated power that is comparable to the star luminosity, but 
this high shear can be maintained for a time of a few years. 

The value of 0.27-0.30 radd -1 we derive for KIC 8429280 
is the largest ever measured, to our knowledge, for a young star 
less massive than the Sun. The analysis of Kepler photometry ac- 
quired in subsequent epochs will enable us to detect any eventual 
variation in the surface differential rotation. 



6. Conclusion 

We have performed an accurate analysis of high-resolution spec- 
tra and high-precision Kepler photometry of KIC 8429280 de- 
signed to characterize the chromospheric and photospheric ac- 
tivity. 

Thanks to the high-resolution spectra we have derived, for 
the first time, the astrophysical parameters, T e s, logg, [Fe/H], 
rotational and heliocentric radial velocity, and lithium abun- 
dance. The lithium abundance allows us to estimate an age of 
about 50 Myr, that, with a mass of about 0.9 Mq, is likely to 
imply a pre-main sequence stage with the star approaching the 
ZAMS along its radiative track. 

A high chromospheric activity level was inferred from differ- 
ent spectral diagnostics (Ha, H/?, Ca 11 H&K, and Ca 11 IRT lines) 
by applying the spectral subtraction technique. Both the Balmer 
decrement, measured as the ratio of Ha to H/3 emission flux, 
and the flux ratio of two Can IRT lines, which is indicative of 
optically thick emission similar to that observed from extended 
facular regions. 

The application of a robust spot model, based on a Bayesian 
approach and a Markov chain Monte Carlo method, to the high- 
precision Kepler photometry spanning 138 days, has allowed us 
to map the photospheric spots and study their time evolution. 
Seven spots were needed to perform a reasonable fit to the Kepler 
light curve, although, because of the exceptional precision of the 
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Kepler photometry, it is impossible to reach the noise floor with- 
out increasing significantly the degrees of freedom and, conse- 
quently, the non-uniqueness of the solution. The active regions 
turn out to be mainly located around three latitude belts, i.e. 
around the star's equator and around ±(50°-60°), with the high- 
latitude spots rotating slower than the low-latitude ones. 

The equator-to-pole differential rotation dQ 0.27 radd 1 
is one of the largest values ever measured for a star cooler than 
the Sun. This is at variance with some recent theoretical models 
that also predict a moderate solar-type d ifferential rotation for 
very rapidly rotating main-sequence stars dKiiker et alj201 lb but 
could fit the scenario proposed by other mode lers of a higher dQ , 
which can c hange along the activity cycle (ICovas et alJl 2005: 
lLanzall2006l) . 
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